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ABSTRACT 

In this paper we show that the most luminous supernova discovered very recently, ASASSN-151h, 
could have been powered by a newborn ultra-strongly-magnetized pulsar, which initially rotates near 
the Kepler limit. We find that if this pulsar is a neutron star, its rotational energy could be quickly lost 
as a result of gravitational-radiation-driven r-mode instability; if it is a strange quark star, however, 
this instability is highly suppressed due to a large bulk viscosity associated with the nonleptonic weak 
interaction among quarks and thus most of its rotational energy could be extracted to drive ASASSN- 
151h. Therefore, we conclude that such an ultra-energetic supernova provides a possible signature for 
the birth of a strange quark star. 

Subject headings: dense matter — stars: neutron — stars: rotation — supernovae: general — super¬ 
novae: individual (ASASSN-151h) 


1. INTRODUCTION 

Fast developing sky survey programs have gathered 
some superluminous supernovae (SLSNe) (Quirnby et al. 
2011; Chomiuk et al. 2011; Gal-Yam 2012) whose peak 
magnitudes are < —21 mag (for a recent review 

on observational properties of SLSNe see Gal-Yam 
2012). SLSNe have been divided into types I 
(hydrogen-deficient) and II (hydrogen-rich), most of 
which cannot be explained by the 56 Ni-powered model 
(Inserra et al. 2013; Nicholl et al. 2013, 2014). Cur¬ 
rently, most of SLSNe-II are explained by the ejecta- 
circumstcllar medium interaction model (Chevalier 1982; 
Ginzburg & Balberg 2012), while the most prevail¬ 
ing model explaining SLSNe-I is the magnetar (ultra- 
highly-magnetized pulsar) model (Inserra et al. 2013; 
Kasen & Bildsten 2010; Woosley 2010), which sup¬ 
poses that a nascent magnetar just after the core¬ 
collapse supernova converts its rotational energy to 
the heating energy of supernova (SN) ejecta 1 . The 
magnetar-powered model has been used to explain a 
few SLSNe-I (e.g., Inserra et al. 2013; Nicholl et al. 2013; 
Howell et al. 2013; McCrum et al. 2014; Vreeswijk et al. 
2014; Nicholl et al. 2014; Wang et al. 2015b; Kasen et al. 
2015) as well as luminous SNe Ic with peak magnitudes 
~ —20 mag (Wang et al. 2015c; Greiner et al. 2015). 

Recently, Dong et al. (2015) observed and analyzed a 
very luminous optical transient at redshift z = 0.2326, 
ASASSN-151h, and found that it reached ~ 2.2 x 10 45 

1 Nonrelativistic pulsar-powered supernovae were origi¬ 
nally proposed by Ostriker & Gunn (1971). Ultra-relativistic 
pulsar-powered shocks were put forward to predict temporal 
plateaus of gamma-ray burst afterglows (Dai & Lu 1998a,b; 
Zhang & Meszaros 2001; Dai 2004). This prediction is well 
consistent with the shallow decay phase (where the flux density 
decays as oc t~ a f with slope of cxj ~ 0 — 0.5) of many early 
afterglows discovered by the Swift satellite (Zhang 2007; Yu & Dai 
2007; Dall’Osso et al. 2011). 


ergs s _1 at time ~ 15 days after the peak. Several expla¬ 
nations of this unusual transient except for the SLSN ori¬ 
gin are implausible. First, the authors excluded a grav¬ 
itational lensing event since the observed redshift is so 
low that any significant contribution from gravitational 
lensing is highly impossible. Second, the active galac¬ 
tic nucleus origin is disfavored because the observed flux 
variability and spectral slope of a normal active galac¬ 
tic nucleus are inconsistent with those of ASASSN-151h. 
Blazars and related “jetted” sources are also excluded be¬ 
cause of their featureless, non-thermal, power law spec¬ 
tra (Dong et al. 2015). Third, a nuclear outburst in 
an active Seyfert galaxy or an event of a supermassive 
black hole tidally disrupting a star is hydrogen-rich, while 
ASASSN-151h’s spectrum is hydrogen-deficient. Thus, 
neither a nuclear outburst nor a tidal disruption event 
can explain ASASSN-151h (Dong et al. 2015). In addi¬ 
tion, ASASSN-151h’s spectrum closely resembles another 
SLSN-I, PTFlOcwr/SN 2010gx. Therefore, Dong et al. 
(2015) concluded that ASASSN-151h is a SLSN-I, which 
is much brighter than all of the other SLSNe discovered 
so far. Dong et al. (2015) argued that if this SLSN is 
powered by 56 Ni, the required mass of 56 Ni is > 30 Mq, 
which is enigmatical for any stellar explosion; if the 
magnetar model is employed, the stellar initial rotation 
period could be ~ 1 ms. A similar estimate for the 
the initial period of the magnetar was also obtained by 
Metzger et al. (2015) and Kashiyama et al. (2015). Con¬ 
sidering the effect of gravitational radiation, Dong et al. 
(2015) further argued that ASASSN-151h challenges the 
magnetar model. 

Provided that it is a newborn neutron star (NS), this 
magnetar is indeed spun quickly down to a period > 
5 ms (Madsen 1998) due to a significant rotational energy 
loss caused by the r-mode instability (Andersson 1998; 
Friedman & Morsink 1998). Thus, the magnetar model 
seems to be ruled out by ASASSN-151h. In this paper, 
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we attempt to resolve this difficulty and constrain the 
physical nature of a magnetar assumed to power such a 
SN. We find that the central engine of this SN could be 
a strange quark star (SQS). 

It was recognized more than forty years ago (Bodmer 
1971) that strange quark matter consisting almost en¬ 
tirely of comparable numbers of deconfined up, down 
and strange quarks could be absolutely stable. If this 
is true, then SQSs composed of such matter could in 
principle exist (Witten 1984) (for an early simple dis¬ 
cussion see Itoh 1970). SQSs were suggested to form 
possibly during some astrophysical processes, e.g. SN ex¬ 
plosions, accretion in X-ray binaries, and mergers of two 
NSs (Gentile et al. 1993; Dai et al. 1995; Cheng & Dai 
1996; Dai & Lu 1998b). Typically, SQSs are under¬ 
stood based on the phenomenological MIT bag model 
of strange quark matter (Farhi & Jaffe 1984), but their 
structural and cooling features are similar to those of NSs 
for stellar masses above one solar mass (Haensel et al. 
1986; Alcock et al. 1986). Thus, it is not easy to identify 
a SQS from observations on stellar structures and sur¬ 
face temperatures (for a review see Weber 2005). Fortu¬ 
nately, a quick rotational energy loss due to the gravita¬ 
tional radiation-driven r-mode instability in a newborn 
(very hot) rapidly rotating NS is absent in a SQS coun¬ 
terpart (Madsen 1998). Therefore, only a newborn SQS 
can reach the Keplerian rotation limit. 

This paper is organized as follows. In section 2, we give 
an upper limit of the initial period and corresponding 
magnetic field strength of a magnetar assumed to power 
ASASSN-151h. In section 3, we fit the light curve of this 
SN. In section 4, we analyze the r-mode instability in the 
magnetar and then constrain the stellar nature. Finally, 
we present conclusions and discussions in section 5. 


2. LIMITS ON MAGNETAR PARAMETERS 

We start by discussing the upper limit of the initial 
period and the field strength of a magnetar (with mass 
of M = Mi .4 x 1.4M 0 and radius of R = R§ x 10 6 cm) 
assumed to power ASASSN-151h. The stellar luminos¬ 
ity due to magnetic dipole radiation is (Ostriker & Gunn 
1971; Shapiro & Teukolsky 1983) 

£mag (t) = ~ (l + t / T )2’ W 

where E Iot = (1/2)/(27t/P;) 2 is the initial rotational en¬ 
ergy with / being the moment of inertia and P, being 
the initial rotation period, and r mag is the spin-down 
timescale of the magnetar and can be expressed by 

T rnag = 2n 2 B 2 R e = da y S ’ (2) 

where c is the speed of light, 1 45 = // 10 45 g cm 2 , Pr ms = 
Pi/1 ms, and B = 10 14 Pi 4 G is the surface magnetic field 
strength. In equation (2), it has been assumed that the 
angle between the magnetic axis and rotation axis is 45°. 

According to the Arnett law (Arnett 1979, 1982), 
which indicates that the peak luminosity L p k of a SN 
at peak time t p k is equal to the instantaneous energy 
deposition rate, that is, L mag (t)|i =tpk = L p k, we obtain 

2 = 27r 2 PP 3 _ 27T 2 P 2 P%k 

‘ (3 c 3 ) 1 / 2 ^ 2 3c 3 / 


From equation (3), we can see that the initial period P; 
reaches a maximum at 


B = B rr = 


(3c 3 ) 1 / 2 / 


= 52.1x 10 14 /45PJ M 


2P 3 Z/ p { f P k 


?-3 T - 1 / 2 -!-! 
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where L pk ,44 = L p k/10 44 erg s \ and t pk ,d = ipk/lday. 
The maximum value of P, is 


P — 

1 i,max — 


t / 1 / 2 


V2Lll%l 2 


= 23.9/ 4 1 fp- k 1 / 4 2 4 t- 1 / d 2 ms. (5) 


Based on the V-Band observation, we can see that 
t P k = i P k,obs/(l + z ) ( z = 0.2326) is ~ 30 to 40 days. 
Around this peak time, the two observational luminosi¬ 
ties derived from the optical band are 2 . 85 / 11/97 x 10 45 
ergs s _1 and 2 .II +/72 x 10 45 ergs s _1 , and at ~ 15 days 
after the peak time, the luminosity is 2 . 16 + 9/9 x 10 45 
ergs s -1 (Dong et al. 2015). Therefore, it is reasonable 
to adopt Lpk ~ (2.0 to 4.0) x 10 45 ergs s _1 . 

The realistic stellar moment of inertia for M > Mq is 
estimated by / ~ 0.237MP 2 (1 + 2.84?y + 18.9?7 4 ), where 
77 = GM/(Rc 2 ) = 0.208Mi.4Pg 4 (Lattimer & Schutz 
2005). Letting Mi .4 = 1 and P 6 = 1 and substituting 
these typical values into equations (4) and (5), we get 
P cr ~3x 10 13 G and Pi. max ~ 0.76 ms. Therefore, we 
conclude that ASASSN-151h could have been powered 
by an ultra-strongly-magnetized pulsar rotating with a 
sub-millisecond period. 

The period of a rotating pulsar has a theoretical 
lower limit, viz., the Kepler period Pk- According to 
Haensel et al. (2009), Pk — C'k M^l / 2 Pg' 2 ms, where 
Ck = 0.78 for a NS and C'k = 0.73 for a SQS. Thus, 
it is easy to find that Pk — 0.78 ms for a NS and 
P K — 0.73 ms for a SQS, provided that M 14 = P 6 = 1. 
These values of the Kepler period basically satisfies the 
following physical requirement: Pk < Pi,max- This 
in fact leads to a new constraint on the stellar mass- 
radius relation (for M > Mq), which can be fitted by 
M > (0.497 + 1.043P 6 - O.113P|)M 0 for a NS and 
M > (0.477 + 0.980P 6 - 0.106P|)M e for a SQS. These 
inequalities rule out very soft equations of state for dense 
matter. 

If the mass and radius of a central rotating pulsar are 
given, we can obtain constraints on the initial rotation 
period and magnetic field strength from observations on 
ASASSN-151h. According to equation (3), Fig. 1 shows 
such constraints based on L p k ~ (3.0±1.0) x 10 45 ergss -1 
and f p k ~ (35 ± 5) days. As a reasonable approxima¬ 
tion, therefore, we adopt Pi ~ 0.8 ms, and B ~ B CI ~ 
3 x 10 13 G, which is significantly weaker than the field 
strengths of magnetars used to power the other SLSNe-I, 
where the latter are generally extremely-strongly magne¬ 
tized, B > 10 14 G (Kasen & Bildsten 2010; Inserra et al. 
2013; Nicholl et al. 2013, 2014). 


3. FIT TO ASASSN-15LH 

We consider a widely-adopted semi-analytical model to 
fit the light curve of ASASSN-151h. Based on the Arnett 
law (Arnett 1979, 1982), letting the initial radius of the 
progenitor Po ~t 0 an( i taking into account the gamma- 
ray and X-ray leakage, the luminosity of a SN can be 
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Fig. 1. — Constraining P[ and B from observations on ASASSN- 
151h for NSs and SQSs, assuming Mi .4 = 1 and Rq = 1. The thick 
solid line and shadow region indicate the best and la values of P[, 
respectively. The blue and red dashed lines are the Kepler limits 
of NSs and SQSs, respectively. The green dashed line is plotted by 
the condition of r ma g > T m , which leads to Pi, ms > 0.247Bi3 (see 
the text in section 3). 


given by 





-^mag 

l”ra 


(*') 



dt\6) 


where r 7 (t) = At 2 = ( 3 «; 7 M e j/ 47 n; 2 c ) £ 2 is the op¬ 
tical depth to gamma-rays (Chatzopoulos et al. 2012), 

T m = (2/tMej/^«; S c) ' is the diffusion timescale, k 7 is 
the opacity to gamma-rays, k is the opacity to optical 
photons, M e j is the ejecta mass, v sc is the scale velocity 
of the ejecta (Arnett 1982), and £ ~ 13.8. In fact, v sc is 
approximately equal to the photospheric expansion ve¬ 
locity u p h, which is estimated to be ~ 2.0 x 10 4 km s _1 
for ASASSN-151h (Metzger et al. 2015). 

We fit the observed data of ASASSN-151h using equa¬ 
tion ( 6 ) and adopting Mi .4 = 1 and Rq = 1 in Fig. 
2. The other parameters are: v sc — 2.0 x 10 4 km s _1 , 
k = 0.1 cm 2 g _1 , M e j = 15Mq 2 , B = 3.1 x 10 13 G, 
Pi = 0.77 ms, and k 7 = 0.03 cm 2 g _1 . In this fit, 
the peak luminosity and the rise time are 3.04 x 10 45 
ergs s -1 and ~ 35 days, respectively. The value of r m is 
~ 31 days, slightly less than the rise time. Theoretically, 
it is required that r mag > r m for magnetar-powered SN 
ejecta to produce sufficiently strong photospheric emis¬ 
sion (Wang et al. 2015a), which leads to a new constraint 
on the initial rotation period and magnetic field strength, 
Pi,ms > 0.247Pi3. This constraint has been shown by 
the green dashed line in Fig. 1. The values of B and P 
adopted in Fig. 2 are in agreement with those estimated 
above. 

Replacing the magnetar input function by the 56 Ni 


2 The optical opacity re of the C+O-dominant ejecta is uncertain 
but taken to be in the range from 0.05 to 0.2 cm 2 g —1 in the 
literature. The diffusion timescale r m of a SN depends on re, the 
ejecta mass, and the photospheric velocity. If the photospheric 
velocity is fixed, the diffusion timescale depends on the values of re 
and ejecta mass. When re = 0.1 cm 2 g 1 , the observed light curve 
requires that the ejecta mass is frozen to ~ 15Mq. To maintain the 
same peak luminosity and shape of a light curve, re X M e j = constant 
must be required. Hence, a smaller (larger) re requires a larger 
(smaller) M e j, i.e., the ejecta mass inferred from the light curve 
depends on the value of re adopted. 



Fig. 2.— Fit to ASASSN-151h using the magnetar- and 56 Ni- 
powered models. Data points are taken from Dong et al. (2015), 
where filled dots are bolometric luminosities derived from the 
optical-ultraviolet band and unfilled dots are derived from the 
optical band. The green dashed, orange red solid, red dot-dash, 
and purple dotted lines are the magnetar input curve, magnetar- 
powered light curve with re 7 —> 00 , magnetar-powered light curve 
with re 7 = 0.03 cm 2 g~ 1 , and 5 ®Ni-powered light curve with re 7 = 
0.027 cm 2 g —1 , respectively. The late-time rebrightening is specu¬ 
lated to result possibly from an interaction between the SN ejecta 
and circumstellar medium, which is similar to SLSN-I iPTF13ehe 
(Yan et al. 2015; Wang et al. 2015d). The vertical blue dashed line 
represents the start of this interaction phase. 


decay input function (see equation 4 in Wang et al. 
2015c), we can reproduce the light curve powered by 
the 56 Ni cascade decay, shown in Fig. 2. For the 56 Ni- 
powered SNe Ic, the fiducial value 3 of tt 7 is usually as¬ 
sumed to be 0.027 cm 2 g -1 (e.g., Cappellaro et al. 1997; 
Mazzali et al. 2000; Maeda et al. 2003). The required 
mass of 56 Ni is ~ 26OM 0 , significantly larger than both 
the lower limit (~ 30M 0 ) estimated in Dong et al. (2015) 
and the mass of 56 Ni synthesized by any possible stel¬ 
lar explosion, demonstrating the invalidity of the 56 Ni- 
powered model. We speculate that the rebrightening in 
the late-time light curve (Brown 2015; Dong et al. 2015) 
could be due to an interaction between the SN ejecta 
and circumstellar medium (where the start of this inter¬ 
action phase is marked by the vertical blue dashed line in 
Fig. 2), which is similar to iPTF13ehe (Yan et al. 2015; 

3 For the C+O ejecta (SNe la and Ic), the opacity to gamma-rays 
is assumed to be gray and does rely on the energy of gamma-rays 
E'y and Y e (where Y e ~ 0.5 is the number of electrons per baryon). 
Therefore it is reasonable to assume the same « 7 for all SNe la and 
Ic with different masses. Monte Carlo simulations performed by 
Swartz et al. (1995) showed that the value of /t 7 is (0.06 zb 0.01)17; 
cm 2 g —1 . When E 1 >0.1 MeV, the value of /c 7 depends only on 
the value of Y e , i.e., « 7 ~ 0.05Ye — 0.025 — 0.03 cm 2 g -1 (see Fig. 
1 in Swartz et al. 1995). Because E^ ( 56 Ni—> 56 Co) = 0.148 MeV 
and E'y ( 56 Co—>- 56 Fe) = 0.847 MeV, we have k 7 ~ 0.025 — 0.03 
cm 2 g -1 . Thus, /c 7 ~ 0.027 cm 2 g — 1 was adopted in the literature 
investigating SNe la and Ic. For complete gamma-ray trapping, the 
total 56 Ni mass is the same as in the case of /t 7 = 0.027 cm 2 g —1 , 
because in the case of a larger /c 7 , the gamma-ray leakage only 
influences the post-peak luminosity. In the case of no gamma- 
ray trapping, however, any amount of 56 Ni does not provide an 
optical-IR luminosity for SNe, since no gamma-ray is trapped and 
converted to optical-IR radiation. 
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Wang et al. 2015d; Dong 2015). Thus, fitting the late- 
time bump is beyond the scope of this paper and we no 
longer discuss it here. 

Since Lp^cx.1 PP max t~^ in equation (5), if r m (~ f pk ) is 
varied, A pk should also change. Decreasing r m (~ i pk ) 
would result in an increase of L pk , but the magnetar 
input curve intersects every peak luminosity. For ex¬ 
ample, when T m is halved, the peak luminosity must 
be ~ 6 x 10 45 ergs s _1 , about one magnitude brighter 
than the peak luminosity of ASASSN-151h (~ 3 x 10 45 
ergs s _1 ). Therefore, it is expected to observe SNe 
brighter than this SN in the future. 

4. R-MODE INSTABILITY ANALYSIS 

In a rotating compact fluid star, the r-mode instabil¬ 
ity is excited by gravitational radiation and suppressed 
by viscosities (Andersson 1998; Friedman & Morsink 
1998; Owen et al. 1998). For a NS, the typical 
driving timescale of gravitational radiation is f gw ~ 
4:7M^IRq 4 P^ s s, where P ms = 27r x 10 3 /fl with f 1 being 
the rotation angle velocity, the dissipation timescale due 
to shear viscosity is t sv ~ 6.7 x 10 7 M^ 4 Rg 3 ^ 4 Tg s, and 
fbv — 2.7 x 1Q 11 MiaRq 1 T 9 _6 P^ 1s s due to bulk viscos¬ 
ity, where Tg is the core temperature in units of 10 9 K 
(Owen et al. 1998). 

For a SQS, the timescale of gravitational radia¬ 
tion is t gw ~ 21M^IRq 4 P^ s s, and the dissipa¬ 
tion timescale due to shear viscosity is t sv ~ 7.4 x 
10 7 dg Q 3 1 M 1 ~ 4 ^ 9 Pg 1 ^ 3 Pg^ 3 s, where a Sj o.i is the strong 
coupling constant in units of 0.1 (Madsen 2000). The 
dissipation timescale t b v due to bulk viscosity as a re¬ 
sult of the non-leptonic weak interaction among quarks 
in a SQS is highly dependent on temperature. The im¬ 
portance of this dissipation mechanism was first recog¬ 
nized by Wang & Lu (1984). For T <C 10 9 K, thv — 
7.9M 1 2 4 P- 4 T 9 - 2 p2 s m( 0 4 0 s, and for T » 10 9 K, t hv ~ 
2.3R§T 2 P 4 s to/ 0 4 s (Madsen 2000), where mioo is the 
strange quark mass in units of 100 MeV. 

An additional viscosity is due to the effect of a 
boundary layer under the solid crust of an old NS 
(Bildsten & Ushomirsky 2000). For a newborn (very 
hot) NS, however, this effect may be unimportant, be¬ 
cause the stellar outer layers are fluidic or at most 
in a mixed fluid-solid state within the first ~ 10 7 s 
(Lindblom et al. 2000). In this early stage, any supercon¬ 
ductivity and/or superfluidity phase in the stellar interior 
does not occur and thus its effect can also be neglected. 

We now study the effect of the r-mode instability on 
the spin evolution of newborn NSs and SQSs. Follow¬ 
ing Ho & Lai (2000), Andersson & Kokkotas (2001), and 
Yu et al. (2009a), we parameterize the evolution of spin 
(f2) and r-mode amplitude (a) as 


dil 

3 a 2 Jn 

f 1 + 1 ) 

N 

iV mag 

(7) 

dt 

A+ 

\ tsv ^bv / 

A+MR 2 ' 

and 





da a 

aA- l 

'! + 1 Y 

aN mag 

(8) 

dt tgw 

A + l 

^tsv ^bv / 

A+MR 2 n' 

where A± = 

/ ± 3a 2 J/2 with J 

= 0.261 and 

J = 


1.635 x 10 2 . The torque of magnetic dipole radiation 


can be expressed by N mag = —P 2 P 6 D 3 /(12c 3 ), where 
the inclination angle of the magnetic axis to the rota¬ 
tion axis has been taken to be 45°. As it grows, the 
r-mode amplitude would eventually saturate, at which 
(da/dt) | Q=Qs = 0 with a s being the saturation ampli¬ 
tude. After then, a s keeps nearly constant and the spin 
evolution follows (Andersson & Kokkotas 2001) 

dfl _ 3 a 2 s Jn 1 N mag 

dt A-(a s ) fgw MR 2 A-(a s ) 

The cooling processes in the magnetar must be known 
to solve equations (7) and (8). As usual, we consider the 
direct URCA processes to describe the cooling of a SQS 
and the modified URCA processes in a NS as well as the 
viscous heating within these two classes of stars. 

For a newborn rapidly-rotating NS, the r-mode am¬ 
plitude is usually assumed to grow to a constant a s 
of order unity (Owen et al. 1998). This order is sup¬ 
ported by numerical simulations (Stergioulas & Font 
2001; Lindblom et al. 2001). Recent studies (Yu et al. 
2009b; Alford et al. 2012) show that a s is in the range 
of ~ 0.1 to a few. Therefore, it is reasonable to take 
a s = 0.1 and 1. In addition, the initial amplitude 
«o = 10 -6 . The spin and temperature evolution can 
be calculated and shown in Fig. 3. We see from the up¬ 
per panel of this figure that a newborn NS rotating with 
~ 0.8 ms quickly reaches the window of instability, then 
the r-mode amplitude grows to a s within hundreds of sec¬ 
onds, and finally this NS keeps r-mode unstable until it 
departs from the instability window at period P ~ 5 ms. 
The green lines (and brown lines) with different mag¬ 
netic fields are overlapped, indicating that gravitational 
radiation dominates the spin evolution. 

For a SQS, the instability window splits into a high- 
temperature part and a low-temperature part due to a 
large bulk viscosity associated with the nonleptonic weak 
interaction among quarks (Madsen 2000). We see from 
the lower panel of Fig. 3 that a newborn SQS with 
Pi ~ 0.8 ms initially lies in the high-temperature insta¬ 
bility region, but this instability window does not affect 
significantly the spin evolution, because the star cools 
so fast that the r-mode instability has no enough time 
(~ 90s) to spin it down obviously. Therefore, when 
it exits from the high-temperature instability window, 
the star has a period close to the initial one. Subse¬ 
quently, stars with different magnetic fields have dif¬ 
ferent spin evolution curves. First, a star with field of 
~ 3 x 10 12 G (thick dashed line), when it cools to a few 
10 8 K, reaches the low-temperature instability window, 
at which time (~ 10 6 s) the star almost keeps its initial 
period, and then it starts to be spun down significantly 
due to the r-mode instability. About 10 8 s later after the 
birth, the star departs from the low-temperature insta¬ 
bility window and loses its angular momentum via mag¬ 
netic dipole radiation. Thus, such a star cannot power 
ASASSN-151h because of a low L mag and r-mode insta¬ 
bility occurring ~ 10 6 s later after the birth. Second, 
after it leaves the high-temperature instability window, 
a star with field of ~ 3 x 10 13 G (thick solid line) is spun 
down quickly by magnetic dipole radiation, so that it 
does never meet the low-temperature instability window. 
Such a SQS, if its initial period P, ~ 0.8 ms, could just 
power ASASSN-151h. Third, a star with much stronger 
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Fig. 3.— Rotation frequency versus temperature for the r-mode 
instability: NSs (upper panel) and SQSs (lower panel) for M 1.4 = 1 
and Rq = 1. The blue, purple and red lines are critical curves 
between r-mode stability and instability regions. The thick dashed, 
solid, and dotted lines show magnetic fields of 3x 10 12 G, 3x 10 13 G, 
and 3 x 10 14 G, respectively. For NSs, the green lines correspond to 
Oi s — 1 and the brown lines to a s =0.1. For SQSs, the purple line 
corresponds to mioo = 2, and the red line to mioo = 1, where the 
thick dashed line gives a s — 0.01 after ~ 10 6 s. The numbers near 
black dots are times since the birth of a SQS. The arrows represent 
evolution directions. 



mode instability. The green and brown lines are plotted for a s = 1 
and 0.1, respectively. All of the lines have considered the effect of 
magnetic dipole radiation. 

field (e.g., ~ 3 x 10 14 G, indicated by the thick dotted 
line) loses its rotational energy so quickly via magnetic 
dipole radiation that it can only drive a SN (with a simi¬ 
lar peak time) much fainter than ASASSN-151h, as shown 
in equation (3). 

We also calculate the contribution of the gravitational 
radiation-driven r-mode instability to spin evolution of 
a NS in detail, which is shown in Fig. 4. We can see 
that for a NS, gravitational radiation will spin down the 
magnetar significantly and bring away most of the stellar 
rotational energy if a s is in the range of 0.1 to unity. 


This implies that a nascent NS seems to be unable to 
drive ASASSN-151h. 

5. CONCLUSIONS AND DISCUSSIONS 

ASASSN-151h is the most luminous SN discovered so 
far (Dong et al. 2015). We have shown that if this SN 
was powered by the 56 Ni cascade decay, ~ 26OM 0 of 56 Ni 
must be synthesized, far exceeding the 56 Ni yield pro¬ 
duced by any single explosion of a massive star. There¬ 
fore, the 56 Ni-powered model is excluded. 

We alternatively adopted the magnetar-powered model 
to explain the high luminosity of this SN. We found that 
the initial period of the magnetar powering ASASSN- 
151h is ~ -fi,max ~ 0.8 ms, which is very close to the Ke¬ 
pler limit. The magnetar’s field is only B ~ 3 x 10 13 G. 
Detailed calculations indicate that a newborn NS with 
these stellar parameters is subject to the r-mode insta¬ 
bility and loses most of the rotational energy via r-mode 
gravitational radiation, but a nascent SQS counterpart 
does not suffer from this instability. Hence, ASASSN- 
151h could have been powered by a SQS rather than a 
NS. In addition, we showed that SQSs with similar pe¬ 
riods but with much weaker or much stronger magnetic 
fields are implausible to explain this SN. 

It should be noted that many other SLSNe-I have been 
assumed to be powered by magnetars with periods less 
than 5 ms, e.g. SN 2011ke (Pi ~ 1.7 ms), SN 2011kf 
(Pi ~ 2.0 ms), SN 2010gx (Pi ~ 2.0 ms) (Inserra et al. 
2013), PTF 12dam (P; ~ 2.6 ms) (Nicholl et al. 2013), 
SN 2013dg (Pi ~ 2.5 ms), LSQ12dlf (P ~ 1.9 ms), 
SSS120810 (Pi ~ 1.2 ms) (Nicholl et al. 2014). If these 
magnetars are NSs, they could quickly lose their rota¬ 
tional energy via r-mode gravitational radiation. Thus, 
while the NS-magnetar model is still challenged in ex¬ 
plaining these SLSNe-I, SQSs as their central engines 
might be able to resolve this question. Therefore, an ac¬ 
curate estimate of the initial rotation periods of newborn 
magnetars powering SLSNe-I would help us distinguish 
between NSs and SQSs, and could eventually provide 
new insights into the stellar nature. 

Although the initial period of the magnetar powering 
ASASSN-151h is approximately equal to the Kepler limit, 
the theoretical upper limit of a SN peak luminosity has 
not yet been reached. We would expect that future sur¬ 
vey programs are promising to discover SLSNe luminous 
than ASASSN-151h 4 . 

Finally, what we would point out is that soft equations 
to state for some exotic matter such as kaon condensation 
in NSs seem to be unlikely, because the maximum mass 
of stars containing such matter is significantly lower than 
the mass measurements (M pu i sar ~ 2.0M 0 ) of two pul¬ 
sars (Demorest et al. 2010; Antoniadis et al. 2013). Even 
if kaons occur in the core region of a neutron star, their 
effect on the r-modes could be insignificant. This is be¬ 
cause the outer regions where the r-modes are mainly 
located do not include kaons, so that we needn’t con- 

4 Dong et al. (2015) estimated that the rate of ASASSN-151h- 
like events is ~ 0.28 — 3.7 Gpc —3 yr —1 (at the 90% confidence 
level), significantly lower than the rate of SLSNe-I (i.e., ~ 11 — 
152 Gpc -3 yr _1 ). Thus, ASASSN-151h-like events are very rare. 
Perhaps SLSNe that are more luminous than ASASSN-151h have 
been missed in previous SN surveys. The cadence of future survey 
programs might be higher than that of current survey programs and 
therefore might detect SNe more luminous than ASASSN-1511 i. 
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sider a bulk viscosity associated with kaons. On the 
other hand, SQSs discussed here are not ruled out if the 
equation of state for strange quark matter is stiff enough 
(Alford et al. 2007). 
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